Observation of a threshold enhancement in the pA invariant mass spectrum 
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An enhancement near the m p +MA mass threshold is observed in the combined pA and pA invariant 
mass spectrum from J/tp — > pK~ A+c.c. decays. It can be fit with an S-wave Breit-Wigner resonance 
with a mass m = 2075 ± 12 (stat) ± 5 (syst) MeV and a width of T = 90 ± 35 (stat) ± 9 (syst) MeV; 
it can also be fit with a P-wave Breit-Wigner resonance. Evidence for a similar enhancement is 
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also observed in ip' — > pK A + c.c. decays. The analysis is based on samples of 5.8 x 10 7 J/ip and 
1.4 x 10 7 ?/>' decays accumulated in the BES II detector at the Beijing Electron-Positron Collider. 

PACS numbers: 12.39.Mk, 13.75.Ev, 12.40.Yx, 13.20.Gd 



An anomalous enhancement near the mass threshold in 
the pp invariant mass spectrum was observed by the BES 
II experiment in J/ip — > jpp decays lj. This enhance- 
ment can be fit with an S-wave Breit-Wigner resonance 
function with a mass around 1860 MeV and a width 
r < 30 MeV .and has been interpreted as a possible bary- 
onium state [2| ■ Similar pp mass-threshold enhancements 
have been observed in the decays B + —y K + pp and B° — » 
D°pp by the Belle Collaboration 0, 0] . These somewhat 
surprising experimental observations have stimulated a 
number of theoretical speculations 0, 0|- It is, there- 
fore, of special interest to search for possible resonant 
structures in other baryon-antibaryon final states. The 
Belle Collaboration recently observed a near-threshold 
enhancement in the pA mass spectrum from B — > pAir 
decays 6]. In this letter, we report the observation of 
an enhancement near threshold in the pA invariant mass 
spectrum in J/ip — > pK~ A and in ip — > pK~ A decays. 
(In this letter the inclusion of charge conjugate modes is 
always implied). The results are based on an analysis of 
5.8 x 10 7 J/ip and 1.4 x 10 7 ip' decays detected in the 
upgraded Beijing Spectrometer (BESII) at the Beijing 
Electron-Positron Collider (BEPC). 

BESII is a large solid-angle magnetic spectrometer 
that is described in detail in Rei.[7|. Charged particle 
momenta are determined with a resolution of u p /p — 

1.78%^/l + p 2 (GeV 2 ) in a 40-layer cylindrical main drift 
chamber (MDC). Particle identification is accomplished 
by specific ionization (dE/dx) measurements in the MDC 
and time-of-flight (TOF) measurements in a barrel-like 
array of 48 scintillation counters. The dE/dx resolu- 
tion is <JdE/dx — 8.0%; the TOF resolution is measured 
to be cttof = 180 ps for Bhabha events. Outside of 
the time-of-flight counters is a 12-radiation-length bar- 
rel shower counter (BSC) comprised of gas tubes inter- 
leaved with lead sheets. The BSC measures the energies 
and directions of photons with resolutions of <je/E ~ 
21%/ v /£(GeV), cj ^ = 7.9 mrad, and a z = 2.3 cm. The 
iron flux return of the magnet is instrumented with three 
double layers of counters that are used to identify muons. 
In this analysis, a GEANT3-based Monte Carlo (MC) 
package with detailed consideration of the detector per- 
formance (such as dead electronic channels) is used. The 
consistency between data and MC has been carefully 
checked in many high purity physics channels, and the 
agreement is reasonable. 

The J/ip — y pK~ A candidate events are required to 
have four charged tracks, each of which is well fitted to 
a helix within the polar angle region |cos(9| < 0.8 and 
with a transverse momentum larger than 50 MeV. The 



total charge of the four tracks is required to be zero. For 
each track, the TOF and dE/dx information are com- 
bined to form particle identification confidence levels for 
the 7T, K and p hypotheses; the particle type of a track 
is assigned to be that of the hypothesis with the largest 
confidence level. In this analysis, reliable identification 
of the K~ is important. To have high efficiency, it is only 
required that one track be positively identified as a pro- 
ton or antiproton. Events where the p, K~ , p and n + 
tracks are all unambiguously identified are subjected to a 
four-constraint (4C) kinematic fit with the corresponding 
mass assignments for each track. For events with ambigu- 
ous particle identification, all possible 4C combinations 
are formed, and the combination with the smallest x 2 is 
chosen. The final \ 2 is required to be less than 20. Fur- 
ther, the p and K~ tracks are required to originate near 
the interaction point, and the invariant mass of the pm + 
combination is required to be less than 1.15 GeV. To sup- 
press background events from J/ip — > pK^YP ', we require 
£ = E miss + 1.39M pK miss < 1.69 GeV (see Fig. [□»), 
where E m i SS denotes the difference between the center-of- 
mass energy (3.097 GeV) and the total energy of the four 
charged tracks, and M p k miss denotes the mass recoiling 
against the proton-kaon system. This selection criterion 
is determined by optimizing the signal to background ra- 
tio based on Monte Carlo simulations. A sample of 5421 
J /ip — > pK~A candidates survive the final selection. The 
pn + invariant mass spectrum for these events, where a 
clear A — > pir + signal is evident, is shown in Fig. EJb). 

The J/ip — y pK~ A events are experimentally quite dis- 
tinct: they contain only charged tracks, three of which 
are heavy particles (i.e., p,p and K), and the kinematics 
strongly constrains the event selection and mass assign- 
ments. In order to maintain a high selection efficiency 
and reduce systematic uncertainties, positive identifica- 
tion of only the proton or antiproton is required, and no 
requirement is placed on the A's secondary vertex. The 
clean A signal in the pir + invariant mass spectrum and 
good agreement with the A signal from a MC sample of 
J/ip — y pK~ A events (Fig. IHb)) indicate that the purity 
of the selected events is very high. 

The level of background in the selected event sample 
was determined with two different MC studies. One used 
a specific set of background processes: J/ip — > pK~T,°; 
AA; AA7T ; pp7r + 7r~; ppn + ir~ ir° ; and all produced 

according to branching ratios from Particle Data Group 
(PDG) Tables || . The fraction of these events that sur- 
vive the J/ip —> pK~ A selection criteria corresponds to 
about 18 events in the selected data sample. The second 
study used an inclusive MC sample of 30 million J/ip 
events generated according to the LUND model |9j . This 
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study predicts that there are 56 background events in the 
data sample. These studies indicate that the background 
in the selected event sample is at the 1 ~ 2% level. 
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FIG. 1: (a) The £ distribution (see text). The solid his- 
togram denotes J/t/j — > pKA events, and the shaded his- 
togram J/ip — > pKTfi events, where both histograms are nor- 
malized to 5.8 x 10 7 J/tp events. Events with £ values below 
the arrow are selected, (b) The pir invariant mass distribution 
for selected events; points with error bars denote the data and 
the histogram the MC (normalized to data). 

The pA invariant mass spectrum for the selected events 
is shown in Fig. [2^a), where an enhancement is evident 
near the mass threshold. No corresponding structure is 
seen in a sample of J/tf) — > pK~A MC events gener- 
ated with a uniform phase space distribution. The pK~ A 
Dalitz plot is shown in Fig. EJb). In addition to bands 
for the well established A*(1520) and A*(1690), there is a 
significant N* band near the K~ A mass threshold, and 
a pA mass enhancement, isolated from the A* and N* 
bands, in the right-upper part of the Dalitz plot. 

This enhancement can be fit with an acceptance 
weighted S-wave Breit-Wigner function together 
with a function /ps(<5) describing the phase space contri- 
bution, as shown in Fig.^c), where fps(S) — N(S 1 ^ 2 + 
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m 



pA 



m p — m-fc, and the param- 
eters a\ and a-i are determined from a fit to the pK~ A 
MC sample events generated with a uniform phase-space 
distribution 
M- 

m = 2075 ± 12 MeV and a width T = 90 ± 35 MeV. 
The fit confidence level is 22.5%( x 2 /d.o.f. = 31.1/26), 
and —2lnL = 29.9. The no resonance hypothesis is also 



The fit is confined to the M p -^ 



M p - 



A < 150 MeV mass region and gives a peak mass of 



tested, and the fit is much poorer: the confidence level is 
5.5 x 10" 10 ( X 2 /d.o.f. = 101.5/29), and -2lnL = 96.2. 
This indicates that the enhancement deviates from the 
shape of the phase space contribution with a statistical 
significance of about la. 

The fit yields N res — 238 ± 57 signal events, corre- 
sponding to a branching ratio 



BR(J/tp 



N res /{2eBR{A 



R-X)BR(X 
pir)) _ 



N 



(5.9 ± 1.4) x 10~ 5 , 



where BR(A -> pn) = 63.9±0.5% is taken from the PDG, 
Nj/^ =L_[5-77 ± 0.27) x 10 7 is the total number of J/ip 
events and e = 5.47 ± 0.05% is the MC-determined 
signal acceptance. 
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FIG. 2: (a) The points with error bars indicate the measured 
pA mass spectrum; the shaded histogram indicates phase 
space MC events (arbitrary normalization), (b) The Dalitz 
plot for the selected event sample, (c) A fit (solid line) to 
the data. The dotted curve indicates the Breit-Wigner sig- 
nal and the dashed curve the phase space 'background', (d) 
The cos 9 P distribution under the enhancement, the points are 
data and the histogram is the MC (normalized to data) 

The signal acceptance and the phase space shape 
fps {$) are corrected for differences between the low mo- 
mentum p and p tracking efficiencies for MC and data. 
The p and p tracking efficiencies are measured with the 
data using a sample of J/ifi — > ppiT + ir~ events. 

A P-wave Breit-Wigner signal function (angular mo- 
mentum L = 1) also gives an adequate fit to the data; 
here x 2 /d.o.f. = 32.5/26 with a mass M = 2044 ± 
17 MeV and a width T = 20 ± 45 MeV, which is con- 
sistent with zero. Fits with higher angular momentum 
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hypotheses L > 2 fail; such states are expected to be 
strongly suppressed near threshold. 

The low acceptance for low momentum protons and 
anti-protons produces a non-uniform acceptance across 
the M pT = 2075 MeV band in the Dalitz plot (Fig.EJb)). 
This is reflected in the non-uniform cos 6 p distribution, 
where 9 P is the decay angle of p in the pA CM frame, for 

-M T < 



TABLE I: Systematic Errors 



the events in the enhancement region [M — Mj 



pA" 

150 MeV), as shown in Fig. Eld) . The distribution agrees 
well with that of a MC sample of J/tp -> KX -> KpA 
with M x = 2075 MeV and T x = 90 MeV. Since the MC 
cos 9 P distribution is generated as a uniform S-wave dis- 
tribution, but the detected MC distribution agrees with 
data in Fig. Hfd) , the observed distribution for the en- 
hancement is consistent with S-wave decays to pA. 

Evidence of a similar enhancement is observed in ip — > 
pK~A, shown in Fig. 01 (a), when the same analysis is 
performed on the ip' data sample. A fit is applied on the 
ip' data sample with the AT(2075) parameters fixed at the 
values obtained from the J/tp data, i.e., Mx = 2075 MeV 
and T x = 90 MeV. The fit shows that the threshold 
enhancement in tp' data deviates from the shape of the 
phase space contribution with a statistical significance of 
about 4.0cr, where the significance is estimated from a 
comparison of log-likelihood values of the fits with and 
without the X(2075) signal function. 




I 6 



M pA (GeV) 



r "•a "! ^ L ". . 

™ ™ H ™ B b b B ^ B 

o □ 

g B B bb °| ■ □ 

OB B °D " BB B 

B □ B b n bm/Vd 00 "b 

[ D B^t° ^^dB^g^B ^ ^ 

4 6 

M^ 2 (GeV 2 ) 



FIG. 3: Results for ip' -> pK~A events: (a) A fit (solid line) 
to the data sample (histogram); the dashed line indicates the 
phase space 'background' contribution. (b)The Dalitz plot. 

The possibility that the enhancement in the J/tp data 
sample is due to interference between N*'s and A*'s has 
been investigated with a partial wave analysis (PWA). 
The PWA results show that if the enhancement were 
from a pure interference effect, many large branching 
ratio J/tp decays to 7V*'s and A*'s near the kinematic 
threshold are required, along with large mutual destruc- 
tive interferences that cancel these large production rates 
[l2|. Also, the similar enhancements seen in the ip' data 
sample and in B — + pA7r observed by the Belle exper- 
iment cannot be due to N* and A* interference effects 
since in these cases, contributions of the signal are far 
from the V* and A* bands in the Dalitz plot (See Fig. |3| 
(b)). 
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Systematic uncertainties from different sources are 
studied. In the above fit, the phase space contribu- 
tion is treated as the 'background' under the enhance- 
ment. Alternative 'background' shape parameters, in- 
cluding iV*'s and A*'s contribution obtained from PWA 
fits, are used to estimate systematic uncertainties from 
the 'background' shape. The fitting bias near threshold is 
checked by MC studies. A set of MC samples combining a 
signal (resonance near threshold) process with a uniform 
phase space process are generated. In each MC sample, 
the mass, width and number of signal events are obtained 
from a fit using the same procedure as that done on the 
data. The averaged offsets between the fit output values 
and their input values are taken as one source of sys- 
tematic uncertainty (fitting bias). The systematic uncer- 
tainty from the tracking efficiencies, especially from the 
low momentum p and p tracks, are checked from data and 
MC comparisons, where the tracking efficiencies of p and 
p are determined from a data sample of J/tp — > ppir + ir~ 
events, and the tracking efficiencies of charged pions is 
obtained from J/ip — > AA, pn events. The systematic 
uncertainty from the kinematic fit is estimated by using 
a different MDC wire resolution simulation model. Sys- 
tematic uncertainties from other sources (such as mass 
resolution) are negligible. The systematic uncertainties 
determined from the above studies are listed in Table [fl 
and the total systematic errors on the mass, width and 
branching ratios are 5 MeV, 9 MeV and 33.5% respec- 
tively. 

In summary, an anomalous enhancement near thresh- 
old is observed in the invariant mass spectrum of pA in 
the J/tp — > pK~ A and ip' — > pK~ A processes. Both S- 
wave and P-wave Breit-Wigner resonance functions can 
fit the enhancement. If it is fitted with an S-wave Brcit- 
Wigner resonance function, the mass is m — 2075 ± 12 ± 
5 MeV, the width is T = 90 ± 35 ± 9 MeV, and_ the 
branching ratio is BR(J/tp -> K~X)BR{X -> pA) = 
(5.9 ± f .4 ± 2.0) x 10~ 5 , where the first errors are sta- 
tistical and the second are systematic. To understand 
the nature of this anomalous enhancement, searching for 
the same enhancement in Kix and i^7r7r modes in the 
J/tp,tp' — > KKtt, KKtttt processes would help to distin- 
guish whether it is from a conventional K* meson or from 
a possible multi-quark state. 
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